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SUMMARY 


Experimental instantaneous pressure records of simulated ram-jet 
models without heat addition showed that they "buzzed" in a manner gov- 
erned internally by quasi -one-dimensional (plane-wave) theory and not 
in accord with the laws of the Helmholtz theory. A theory applicable 
to the buzzing problem was obtained by modification of the one -dimensional, 
unsteady-flow theory by introducing specified assumptions which greatly 
facilitated the computational efforts without impairing the usefulness 
of the theory for this flow. Because of the difficulty in computing the 
three-dimensional unsteady viscous flow when the shock was in motion 
ahead of the cowling and in computing when the shock was in motion far 
from equilibrium with the waves moving upstream from the nozzle, it was 
necessary to use a few experimental points at the start of the cycle as 
a basis on which to compute the many remaining pressure points of the 
cycle. Theoretical computations of the pressure-time curves, which were 
made on this basis for one supersonic-diffuser configuration with two 
different subsonic-duct lengths, agreed closely with the experimental 
records. The entrance and exit mass-flow — time curves, as well as the 
wave diagrams, were also computed for these configurations by the 
unsteady-flow theory. A possible simulation of the buzzing in ram jets 
with combustion by the test of a cold-flow model in which a contraction 
in area replaces the combustion chamber is also discussed. 


INTRODUCTION 


The buzz or pulsation of axially symmetric supersonic inlets having 
central bodies has been the subject of investigations since Oswatitsch 
first observed the phenomenon in 19*1-4 (ref. l). This buzzing is essen- 
tially an oscillation of the normal shock along the spike of certain 
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diffusers operating at mass flows below design, and this buzz decreases 
the performance of the propulsion units appreciably (ref. 2). When a 
fixed-geometry ram jet is accelerated to the design Mach number, it is 
necessary to fly through a range of mass flows less than the design 
value. Consequently, the importance of either avoiding the conditions 
causing the buzz or reducing and controlling the amplitude of the oscil- 
lation is evident. In order to accomplish this end it is, of course, 
desirable to understand the complete nature of the phenomena. Antonio 
Ferri and Louis M. Nucci (ref. 3) have described one of the mechanisms 
which initiate the buzz. This paper continues the analysis of the buzzing 
cycle by means of a modified one-dimensional approach. The buzzing 
referred to herein is not to be confused with the pulsations of inlets 
sometimes observed due to improper operation of the combustion system. 

The causes of this buzz, in the absence of combustion, have been 
studied previously and include effects due to the presence of separation 
inside the entrance to the cowling. The choking resulting from this 
separation reduces the internal mass flow and forces the normal shock 
to move forward; thus, the buzz is initiated. Reference 3 has shown 
that this separation is often due to a vortex sheet near the inner wall 
of the cowling lip. This vortex sheet arises from the intersection of 
either the conical and normal shocks or, in some cases of separation on 
the central body giving rise to lambda shocks, from the intersection of 
the legs of the lambda shock. In certain instances, the separated flow 
from the central body has appeared sufficient in itself to choke the ^ 

cowling and precipitate the pulsations. As the buzzing cycle progresses, 
the normal shock moves rearward and assumes a position where it may remain 
for an appreciable time, depending on internal-flow considerations, before »> 

the start of the next oscillation. 

Some very general theories based on internal mass-flow consider- 
ations have been expounded which give little insight into the aerodynamic 
nature of the phenomena making up the complete buzz cycle (refs. 3 and. 

4). An attempt has also been made (ref. 5) to predict the frequency 
only of the buzz by considering a ram jet to be similar to a Helmholtz 
resonator. The application of a quasi-one-dimensional theory to the 
buzz problem would seem proper from considerations of the ram-jet geom- 
etry. First of all, the area variations are gradual because of the 
necessity for high diffuser performance. Furthermore, the ratio of 
maximum area of the unit to the minimum area usually varies from 3 to 5j 
whereas for the passages of the internal-flow system the ratio of the 
length to the maximum diameter is generally in the range of 6 to 15 . 

These proportions suggest the use of an unsteady-flow theory, such as a 
quasi-one-dimensional one, which assumes a slight variation of the moving 
pressure waves from the plane form, rather than a Helmholtz type theory 
based on two- or three-dimensional waves (that is, large, abrupt area 
changes) . 
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Consequently, an appropriate modified one -dimensional nons tat ionary - 
flow wave theory has been derived at the Langley laboratory. This theory 
assumes quasi-steady flow for the reflections of upstream-moving waves 
which strike the normal shock and for waves passing through an area 
change. Because of the difficulty in computing the three-dimensional, 
nonsteady phenomenon which occurs as the normal shock moves rapidly out- 
ward and inward during part of the cycle, it is necessary to use a few 
experimentally determined values of the waves generated by this motion 
as a basis on which to compute the remainder of the cycle. For the par- 
ticular cycles computed in this paper, four or five experimentally deter- 
mined instantaneous pressures are used to compute one hundred or more 
additional pressure points required to complete the cycle. Once the 
initial problem has been specified by the experimental points, however, 
the remainder of the solution is easily obtained from graphs with a mini- 
mum of computation and good degree of precision. 


SYMBOLS 


a local speed of sound 

A area measured normal to ram- jet axis 

CpCn characteristic parameters for steady flow 

Cp coefficient of specific heat at constant pressure 

Cy coefficient of specific heat at constant volume 

h state ahead of typical downstream wave 

j state behind a typical downstream wave 

k state behind reflection from a typical downstream wave 

L distance from cowling inlet to exit nozzle of ram jet 

m mass flow 

uioo mass flow at infinity through a stream tube of diameter 

equal to cowling lip diameter 

M Mach number 

p local static pressure 
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characteristic parameters for unsteady flow 
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Subscripts: 

b 

c,d 

e 

f ,6 


characteristic parameter reflected at cowling by parameter 
Qj for case when both parameters are evaluated in fluid 

having entropy equal to "effective value" 
gas constant 

nondimens ional entropy referred to free-stream datum 
time, sec 

local static temperature 

energy added to flow across flame front 

local fluid speed 
coordinate along ram- jet axis 

ratio of specific heats, assumed equal to 1.400 

distance from flame front or area contraction to exit nozzle 

angle included between ram- jet axis and line from apex of 
spike to cowling lip 

local static density 

state in ram jet at start of buzz 

states on opposite side of typical entropy discontinuity 
effective one-dimensional values 

states on opposite side of entropy discontinuity entering 
exit nozzle 

station at cowling inlet (total cross section enclosed by 
inlet) 
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n exit nozzle (measured, not effective, value) 

q state "behind an upstream wave generated by entropy discontinuity 

passing through nozzle 

r state behind second entropy discontinuity entering exit nozzle 

00 free^stream state 

a region between flame front and nozzle of ram jet 

3 region between area contraction and nozzle of ram jet 

0 state obtained by isentropic deceleration from static con- 

ditions to zero velocity 

1 annular station at cowling inlet 

2 station immediately behind cowling in a simplified ram jet 

3 station immediately ahead of nozzle in a simplified ram jet 

Except for the cases noted in symbol list, subscripts denote the 
evaluation at the state or station designated by the subscript. A bar 
over a symbol indicates an average over a cycle. 


THEORY 

Basic Theory 


One-dimensional unsteady-flow characteristic theory is governed by 
the following differential equations, which determine the variation of 
the parameters P and Q, respectively, along the corresponding character- 
istic lines having slopes of — = u + a and — = u - a 

at dt 


SP = __JL_ D lQ ge A + 7 - 1 a D I , 1 a 5 I 
5t a Ooo Dt y a o 00 Dt 7 a Oc0 St 


SQ 

St 


a D logg A y 


Dt 
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O St 

w oo 


(2) 
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( 3 ) 


Q = 


2 a 

y - 1 a °„ 


u 
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where + (u ± a) is the 

6t bt bx 

, D( ) b( ) . b( ) 
characteristic and — = -r— — + u 

Dt ot 

ular fluid particle. 


< 3 x 


total variation along a 
is the variation of a partic- 


These equations are derived in appendix A by combination of the 
equations of energy, momentum, continuity, and state. Equations of 
similar type may be found in reference 6 in which, however, the area is 
not considered to be time dependent. 


The parameters P and Q, have much the same relation to unsteady- 
flow theory as the parameters Cj and Cjj have to the Prandtl- 

Busemann characteristic theory for steady flow (for example, see refer- 
ence 7). In the two-dimensional steady-flow system, Cj and Cjj are 

invariant along characteristic lines of the first and second family 
except when the characteristics pass through rotational flow fields. The 
parameters Cj and Cjj may vary in an irrotational three-dimensional 

flow because of the effects of stream-tube-area change. Similarly, in 
one-dimensional unsteady flow, the values of P and Q are constant 
along the characteristics except when influenced by entropy changes or 
area changes (the latter factor influences the parameters by accounting 
for the variation of the flow from the exact one -dimensional flow) . In 
steady flow the inclination of the streamlines is proportional to the 
differences of Cj and Cjj whereas the velocity is proportional to 

their sum. Again there is a similarity in one -dimensional unsteady flow 
where the particle speed is proportional to the difference of P and Q 
whereas the local speed of sound is proportional to their sum. 

In the use of the one-dimensional characteristic equations above, a 
very lengthy and tedious step-by-step process is required to obtain a 
solution. Consequently, certain simplifications and assumptions are 
made to obtain a modified theory which appreciably facilitates the pro- 
cedure necessary to obtain a solution. 
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Modified Theory 

In the interest of simplicity and computational facilities, the basic 
theory has been modified by application of the following assumptions. 

Substitution of discontinuous area changes for area gradients . - 
In order to simplify the solution, gradual variations in cross section 
are accumulated and treated as discontinuous area changes. Thus, the 
actual ram jet, as shown in figure l(a), is assumed to have the shape of 
figure l(b ) where the only area changes appear at the cowling and the 
nozzle exit. The penalty for a considerable reduction in computation is 
a time distortion in the wave pattern since both the duration and time 
of generation of reflections is greatly reduced. In other words, the 
unmodified solution yields gradual reflections as soon as an area 
gradient is encountered by a wave. No reflections, however, can arise 
in the simplified ram jet until a wave hits either end of the duct, 
whereupon stronger reflections, representing the accumulation of the 
aforementioned gradual reflections, are generated. 

Restriction requiring equality of entropy of a fluid particle at 
both entrance and exit of hypothetical ram jet .- In order to simplify 
the procedure further, all waves moving through the duct proper are 
assumed to be isentropic. Thus, the entropy increase in the fluid is 
assumed to be the cumulative result of the conical shock, normal shock, 
and separation on the central body preceding the main duct. The entropy 
of entering fluid particles varies only as these factors vary and then 
this entropy (which may have a different value for various particles) 
remains constant for a particular fluid particle flowing along the duct. 

Use of quasi-steady flow to obtain reflections at duct exit and 
entrance .- The assumption of quasi-steady flow is introduced to provide 
a simple solution for the integrated form of the characteristic equa- 
tions when waves pass through area changes or encounter the normal 
shock. This modification is based on the fact that steady flow is 
established at an area change in a very short time interval after the 
passage of the initial wave. Consequently, all the minor transmitted 
and reflected waves may be accumulated with the major transmitted and 
reflected waves, and then the group treated as consisting of only one 
transmitted and one reflected wave. Of course, this treatment also 
yields a wave pattern with a slight time distortion, but, in view of 
the labor saved, this distortion is justifiable. 

The reflection condition at the nozzle exit is then easily deter- 
mined. Since the subsequent analysis is to be limited to supersonic 
flight for which the pressure in the duct is always greater than the 
stagnation pressure required to produce sonic velocity by isentropic 
expansion to atmospheric pressure, the boundary condition of sonic 
velocity in the nozzle is dictated. Consequently, since the ratio of 
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the area of the duct to the nozzle area is invariant, there is auto- 
matically imposed the condition of constant Mach number (-which must be 
equal to the effective Mach number at the start of buzz) in the ram jet 
just ahead of the nozzle under the assumption of quasi-steady flow. 
Thus, from the definitions of P and Q, one finds for the values of 
these parameters existing just ahead of the nozzle 


P 3 ~ Q 3 

M e = — = 2 

a e P3 + Q3 

4 

7-1 


(5) 


or 




- M e 

+ Me 


5 - M e 
5 + M 


7 


1.400 


( 6 ) 


where M e is the effective Mach number of the duct. Thus, the value 
of the parameter P^ for a characteristic moving downstream at the 
nozzle uniquely specifies the parameter of the reflected-upstream 

characteristic. 

For the case of reflections from the cowling inlet of upstream- 
traveling waves, the parameter P can be evaluated if the mass flow 
entering the cowling is known. During a large part of some cycles the 
normal shock is inside the cowling and therefore the mass flow may be 
determined from conical-flow theory and the external geometry. For 
that part of certain other cycles when the normal shock is ahead of the 
cowling but has not progressed out as far as the point where buzz starts, 
experimental steady- flow values of the mass flow can be used. Once the 
value of the mass flow is assumed, then there exists a relation at the 
cowling entrance between P, Q, the entropy S/R, and the effective 
area A e for quasi-steady flow. As an illustration, assume that the 
external geometry is such that the conical shock intersects the cowling 
lip (that is, diffuser operating at design M^) while the normal shock 
is swallowed. Then the following equations of continuity, energy, and 
so forth, apply for quasi-steady flow: 


A iPoAo a oo = A eP 2 M 2 a 2 


(7) 
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a„ = a 


°2 


(8a) 


7 - l/ u c 
2 a. 




= 1.0 


( 8 b) 


Introduction of the entropy term into equation (7) leads to 


7+ 1 e e 

S 2 -b c 

^2/^2\ 7 " = e R 
/ A P 

00 \ 00 / t. 


(9) 


In terms of P and Q equations ( 8 b) and (9) become 


7 - ife ~ Q 2 


) + f 21 ^) ( P 2 + Q 2) 2 - 1 


( 10 ) 


(h*)[<'. * ^-Ar 1 ■ ft) 


7+1 
v?" 1 A, 


S 2- S co 


M — e R 
A 


(11) 


Equations (10) and (ll), plus the restriction that u > 0, which 
requires P > Q, uniquely determine the value of P for a given Q 
and S/R, or vice versa. 

Before the reflected P^’ can be determined, however, the value 
of Q 2 must be found from the corresponding value of along the 

same upstream characteristic (fig. 2). This value of Q 2 may be 
obtained from equation (2) which states, for A = constant and for 


isentropic flow of each elemental fluid volume 


VDt 


- 0 , 


Q 3 - Q 2 -J 2 y a, 


1 a d S = 1 


R 7 


\ a AS 
Z_ a 0 R 


( 12 ) 
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For the low Mach numbers encountered in a ram- jet diffuser 
equation (12) becomes 


a 

a Ooo 




1 , 


Similarly, 


Q: 


= 7 a 


a. r 



1 s 3 - s 2 
7 R 



1 3 4 S 3 ~ S 2 
7 R 


(13) 


(14) 


Therefore, since — is known from the starting conditions and, for a 

R 

So 

given value of P 2 and Q 2 are specified, one can determine from 

equations (10), (11), (13), and (l4) the value of the reflected charac- 
teristic parameter py. The detailed procedure to find Py for a 

specific Q3 would be the following (see fig. 2): 

(1) Equation (10) is solved for P 2 in terms of Q 2 . 

(2) This value of P 2 is substituted into equation (ll) so that it 

becomes an equation of the form: Q 2 = Q 2 





(3) This equation for Q 2 is then plotted together with equation (13), 

which is also in the form Q 2 = Q 2 ^-^j, an d the intersection of the two 
curves determines the correct value of S 2 /R and Q 2 . 

(4) The value of P 2 is then found from Q 2 and equation (10). 

(5) The term Pj' is then evaluated by equation (l4). 

The most convenient method of obtaining the values of P^ ' from 
Qj, however, has been found to use a plot of Py against where 

py and Q5 are computed by adjusting each point of the curve for 

against of equations (10) and (ll) by a factor of 

1 S 3~ S 2 

7 R 
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This manner of determining the reflection at the forward section 
of the ram jet is inapplicable when the motion of the normal shock is 
far from a quasi-steady, or near equilibrium, state with the waves which 
have come upstream from the nozzle. In other words, when the shock has 
been set in motion from a point on the external spike owing to three- 
dimensional waves emanating from the cowling and then the shock enters 
the cowling, it may not be in equilibrium with the one-dimensional waves 
existing inside. Another case where the method is inapplicable would be 
the time when separation generated at the foot of the shock causes the 
shock to move, regardless of the upstream waves from the throat. 

In addition, the use of some three-dimensional unsteady- flow char- 
acteristic theory in the presence of viscosity would be necessary to 
determine the reflection pattern from the normal shock when it is in 
motion ahead of the cowling. For these reasons, experimental pressures 
are used to prescribe the downstream waves generated when the normal 
shock is in motion either ahead of the cowling or else too rapidly for 
quasi-steady flow. In the cycles computed in this paper, both condi- 
tions occur simultaneously and necessitate the specification of only a 
few waves. 

Reference of characteristic parameters to the basic entropy and 
subsequent neglect in wave diagram of entropy discontinuities in duct .- 
Although the variation with time of the entropy of the fluid entering 
the duct is allowed, the sweeping modification is employed of referring 
all the characteristic parameters to the basic effective entropy (that 
is, effective entropy of the fluid present in the duct at start of 
buzz). This assumption drastically reduces the computational effort, 
since it eliminates the determination of an entropy-discontinuity 
streamline every time a wave strikes the normal shock and, then in turn, 
the determination of reflected and transmitted waves each time a wave 
impinges on these entropy discontinuities. It is also consistent with 
the assumption that % 1 in equations (13) and (l4) relating the 

a Ooo 

increment of P or Q to the entropy variation across a discontinuity, 
since that assumption also implies the neglect of waves reflected from 
such discontinuities. 

The use of a basic entropy simply amounts to the insertion of an 
imaginary lamina of fluid having the "basic entropy" in a fluid of 
another value of entropy in order to determine the pressure and velocity, 
but not the density, of the flow. This procedure is permissible because 
the last terms in equations (l) and (2) represent the increment in the 
characteristic parameters P and Q across an entropy discontinuity 
when the pressure and velocity are constant. This fact may be proven 
as follows with reference to figure 3 * By applying the definitions 
of P and Q to the characteristics in figure 3 > the following equa- 
tions are determined: 
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since u c = u^. 
Therefore, 


AP = Aft = 
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P 
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a O 00 
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Introduction of entropy into equation (17) yields 


AP = Aft = 



7-1 AS 

7-1 


2 a c 

e^ R ( 
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2 a o 
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i/V - l csf* . . . 

7 -1 a o M 

27 R 

2\ 27 

R ) 

IAS for 

a c ^ -1 
x 

, 1 


7 R 

a o 

7 R 


in figure 3(h) 'the 

true values of Q and 


(18) 


since p d = p c (19) 


( 20 ) 


entropy were determined as specified by equations (13) and (l4) 

R 


fa 


from the corresponding values in the fluid of entropy -g“- by adding 
S _ go 

the correction R . Consequently, when they are reduced to the 

So 

basic entropy — of the imaginary lamina inserted in the fluid, the 
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factor - 0 


7 R 

and the original values of 


^that is, i ^ from equation (20)^ must be added, 


and Qo are recovered. 


The relationship between the pressure and sonic speed within a 
fluid of constant entropy is: 


P 


P 


e 




( 21 ) 


Since the pressure and velocity, but not density, are invariant 
across an entropy discontinuity in the flow, it is permissible to ignore 
the variation in entropy and use the values of P and Q based on the 


S3 

effective entropy 

t\ 


at the start of buzz for the computation of pres- 


sure as 


27 

— = / R + Q V and velocity as 

p e V P e + <W 


P - Q 
2 


in any part of 


the fluid regardless of entropy. The wave diagram will then be greatly 
simplified, but will not be a true picture of the flow in regard to 
stream density or mass flow. 


Of course, if the entropy at a point is known, the mass flow can 
be determined from the values of P and Q in the imaginary lamina of 
basic entropy. When reference is made to figure 3(b) where condi- 
tions c correspond to the imaginary lamina. 


m c _ Pc u c 
m d ' Pd u d 


1 7-1 S d“ S c 



S d" S c 7-1 


= e 


since p = p 


u = u. 


( 22 ) 


(23) 


(24) 
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When equation (24) is expanded and higher order terms are 
neglected, 


s c _s d 7-1 



Then equation (26) becomes 


^ = Ac * QcV 'V pc - a c \ 

“b \ p b + %1 \ P b " %) 


( 28 ) 


Thus the ratio of tne actual mass flow to the mass flow at the start of 
buzz becomes. 


m d 

si 

1 

m 

c 

m b 

m c 

m b 


/ P c + Q cV V p c ~ Q c \A ^ 7-1 S e ~ 
l P b + %) \ P b - %){ 7 R 


(29) 


There still remains the problem of determining the local entropy. This 
task, if it is to be precisely accomplished, requires the plotting of 
the entropy-discontinuity space-time histories. Without these histories, 
however, it is possible to obtain easily a so-called extrapolated value 
of the mass flow entering the cowling by assuming quasi-steady flow from 
infinity to the cowling lip. This value is not a true value when the 
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shock is inside the cowling and out of equilibrium with the upstream 
Q waves, since the mass flow entering must be equal to m M when the 
conical shock is at the lip of the cowling. This error is due to the 
fact that the parameters P and Q were determined experimentally at 
stations to the rear of the shock and included the effect of the shock 
motion. Consequently, use of these P and Q terms evaluated 

at — = 0 give the mass flow just behind the shock instead of at the 
L 

cowling entrance so that the values are only correct when the shock is 
ahead of the cowling or when it is in equilibrium inside the cowling. 
Since these parameters determine the local pressure and velocity, the 
assumption of quasi-steady flow permits the determination of the local 
temperature from the energy equation. Then, when the pressure and the 
temperature are known, the local entropy and, in turn, the mass flow 
may be found. 

When the presence of entropy discontinuities in the flow is dis- 
regarded in the wave diagram, the waves originating from the arrival of 
these discontinuities at changes in cross section are, of course, 
omitted. In particular, the upstream wave generated by the passage of 
a discontinuity through the choked exit nozzle is neglected and this 
omission must be justified. Figure 4 shows the wave pattern in the 
diffuser created by the passage of an entropy change through the nozzle. 
The following relations then apply 



(30) 
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(34) 


(35) 


(36) 


Thus, waves of strength Qq - Qg are omitted in a solution neg- 
lecting entropy discontinuities ^ which is equivalent to the assumption 

w (fg^M«iV In cases where this assumption is not justified, 
e R / 

s — s 

hut another entropy discontinuity r ~ having opposite sign so that 


^r _ ^f « 1 S oon follows the initial discontinuity through the area 
R R 

change, then the assumption has only introduced a time distortion in the 
wave pattern by treating a positive wave followed by an approximately 
equal negative wave as a single wave of zero strength. 


■Rffective values employed in the modified theory.- One-dimensional 
fl ow presumes no variation ol r the flow properties in a direction normal 
to the stream velocity. Consequently, in order to apply such a theory 
to an actual ram jet where boundary- layer separation and centrifugal 
force actually cause a variation in the fluid properties across a sec- 
tion of the flow, effective values of the fluid properties and cross- 
sectional area must be employed. The effective values are defined as 
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the values satisfying the one-dimensional channel-flow relationships 
regarding continuity, energy, momentum, and the equation of state. 

Since the effective area and Mach number Ag and Me have such 
an important role in the theory, the equations which determine them are 
presented. Steady flow is assumed to exist for an infinitesimal time 
before buzz is started. 

The mass flow at start of buzz is 



(37) 


and 



(38) 


From the definition of entropy and equation of state. 


2 S e” S «> 


r 



(39) 


Since a, 


a, 


*o ( 


'e 


at start of buzz, equation (39) becomes 



(40) 


Combining equations ( 38 ) and (40) yields 


7+1 



( 41 ) 
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Insertion of energy equation into equation (4l) gives 


7+1 



The specification of A e plus experimental values for the other terms 
on the right-hand side of equation (43) prescribes M e uniquely. 

Approximation for simplifying con struction of characteristic 
diagram.- Instead of the step-by-step procedure of constructing the 
characteristic lines with local slopes dependent on the sum or differ- 
ence of the local fluid speed and sonic speed, the value of these speeds 
at the ends of the duct is used to obtain an average slope. Thus char- 
acteristic lines of the first family are drawn with a slope 

u , a 

d 2 + i“ 

L _ °oq u °° (44) 

dt L 


evaluated as the average of the velocities for = 0 and jr = 1. 
Similarly, for the second family. 



( 45 ) 
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whereas for particle space-time histories 



dt _L_ 

a Ooo 




A linear variation of the parameters is used for interpolation between 
characteristics of the same family. 


The characteristic theory has been greatly simplified so that the 
solution is easily obtained with the help of three plots, namely, 


(a) pressure-sonic velocity, 


2 

p - Pe^r 1 , 0>) 


P - Q reflections at 


cowling inlet, and (c) Q - P reflections at the sonic nozzle. (See, 
for example, appendix B and figs. 25 and 26 for typical solutions using 
these plots.) After the experimental values of the few waves generated 
by the external motion of the normal shock are plotted, the remainder 
of the solution consists of reading values off the curves and plotting 
the characteristics. 


APPARATUS 


The apparatus employed for the experiments is shown in figure 5 
and is essentially the same as that used in reference 3 with certain 
modifications necessary for the conversion from steady- flow to unsteady- 
flow measurements. The total-pressure rake at the survey plane 
(x = 11.50 inches), figure 5> was removed and an instantaneous static- 
pressure pickup employed in its place. The instantaneous pressure 
pickups were NACA miniature electrical pressure gages mounted flush 
with the diffuser wall. Also, since only reflections from the nozzle 
exit were desired, the mass-flow orifice was eliminated. 

Tests were made at a Mach number of 1.9 in the Langley 4 -inch by 
5- inch blowdown open jet by using low-humidity air from pressurized 
tanks. The test Reynolds number based on cowling- inlet diameter was in 

6 6 

the range between 2.9 x 10 and 3.6 X 10 . Three different subsonic- 
duct lengths of 34 inches, 156 inches, and 390 inches were employed with 
instantaneous pressure pickups located along the diffuser as shown 
schematically in figure 6. The 90 ° elbows necessary to attain the 
desired duct lengths within the space limitations of the laboratory are 
indicated. The length of the central bodies (which had the same after- 
section in all cases) relative to diffuser lengths is also shown. 
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Tests were made on the three supersonic-diffuser configurations 
listed in the table below: 


Config- 

uration 

Semi cone 
angle 
(deg) 

Maximum diameter 
of central body 
(in. ) 

C owling 
(reference 3) 

Values of cowling- 
position parameters 

( e z) 

I 

25 

0.80 

E 

44° 0», 40° 40*, 36° 50’ 

II 

25 

1.10 

D 

43° io», 39° 50*, 36° 50* 

III 

30 

0.80 

E 

48° 30*, 43° 40*, 40° 40‘ 


The dimensions of these central bodies and cowlings are given in 
reference 3 with the exception of the central body of configuration I, 
the dimensions of which are given in figure 7* 

A commercial concentrated arc light served as the source for a 
shadowgraph system in which the image on a ground-glass screen was 
recorded on either of two l6-millimeter motion-picture cameras. For 
the longer diffuser lengths a low-speed camera operating at 64 frames 
per second was used, whereas for the short diffuser a Wollensak Optical 
Company Fas tax camera running approximately 500 frames per second was 
employed. Signals from the pressure gages were amplified by a 
Consolidated type 1-113 amplifier and then recorded on sensitized paper 
in a Consolidated type 5-114 recorder employing galvanometer elements. 
The frequency response of this system was flat from 2 to 300 cycles per 
second. In order to correlate movie film with the pressure records, a 
timing light was attached to one edge of the shadowgraph glass screen 
and the circuit energizing this light connected in parallel with one of 
the recording galvanometers so that every time the light flashed a 
"blip" appeared on the pressure record. A synchronizing device was 
attached to the low- speed camera which flashed the light on the shadow- 
graph screen once each frame, while 60 - cycle alternating current was 
used to operate the light when the Fastax camera was used. In order to 
reduce the labor of correlation, an interrupter was placed in the light 
circuit to stop the flow of current for a noticeable period four or five 
times a second, and hence provide blank spaces on the film and pressure 
record which could be easily counted. 

A manually operated, 2- inch gate valve served as the exit throttle. 
Throttle position was recorded on the pressure film by means of an elec- 
trical circuit employing a ten-turn linear potentiometer attached to an 
extension of the valve stem. 
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TEST PROCEDURE 


Tests of the various configurations were made by the following pro- 
cedure: the desired central-body projection was set by scale measure- 

ment and, with no tunnel air flowing, a trace was then recorded on the 
pressure film; in the cases of runs made with the $ 4 — frames -per-second 
motion-picture camera, a strip of "no flow" motion-picture film was also 
made. The exit throttle was opened to such an extent that the mass flow 
would be limited only by the supersonic-diffuser configuration. Then 
the tunnel was started and brought up to operating conditions. The motion 
picture camera and pressure recorder were started and the throttle closed 
very slowly until buzz began, at which point the throttle position was 
held approximately constant for a few cycles before continuing on to a 
lower valve opening. This procedure of a slow throttling interrupted by 
pauses at constant valve area was continued for the duration of the test 
run except for a rapid opening of the throttle at the end of the run. 

The tank stagnation pressure as well as the differential-gage reference 
pressure were recorded during the run. A no-flow pressure trace was 
also made after a run when possible. The differential gages and valve- 
position recorded were calibrated statically each day on which tests 
were run. 

The cowling-position parameter, or angle between the axis and a 
ray from the cowling lip to the apex of the central body, was determined 
from enlargements of the no- flow motion pictures. The outlines of the 
conical body were extended to intersect at the apex, and then angles 
were measured from this point to the outline of the cowling lip. This 
method, however, has a low accuracy (of the order of 20 to 30 minutes) 
because of the presence of a slight fuzziness on the shadowgraph 
outlines. 

The value of A n , area of the valve opening normal to the axis, 

was found for each position of the throttle wheel by integrating the 
shadowgraphs obtained by placing the valve in the optic system with the 
flow axis alined with the optic axis. The value of A n is only a rela- 
tive, and by no means exact, measure of the effective throttle area 
since the valve seat is appreciably recessed. Consequently, the effec- 
tive area is always much larger than A n , but the variations in A n 

still indicative of variations in the effective throttle area. 


are 
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RESULTS AND DISCUSSION 

Comparison of Theoretical and Experimental Results 


Theoretical computations of buzzing cycles for configuration I 
(e z = 44° + 20*) with duct lengths of 156 inches and 34 inches have been 

made. These cycles were typical of the configurations which buzzed with 
the vortex at the lip and, in addition, were the most complicated in 
the sense that the one-dimensional wave pattern included a large number 
of generated and reflected waves. The cycle for the 156 -inch duct is 
discussed first because it more clearly shows the waves and their 
reflections since the duct is long enough to permit the complete gener- 
ation of a wave before its own reflection from the nozzle exit returns 
to alter the generation. 

The basic computations for this cycle are given in appendix B. The 
characteristic or wave diagram, including one space-time history for a 
fluid particle which enters as the buzz starts, is shown in figure 8 , 
whereas the theoretical pressures are plotted on the experimental pres- 
sure record in figure 9 . Flagged symbols are employed in this paper to 
denote those values computed from the experimental records. A compari- 
son of the many additional theoretically computed points with the experi- 
mental pressure record shows generally good agreement for all five pres- 
sure traces and thus offers strong evidence that the modified theory is 
valid. 

The lowest trace of figure 9 is that of the valve-position indi- 
cator. This is not a smooth curve because of the electrical "pickup 
and, consequently, the value employed to determine A n is the average 

of the maximum and minimum of that part of the trace where the "pickup 
is smallest. The topmost blips indicate the timing-light flashes used 
to synchronize the motion pictures to the pressure record. 

A strip of movie film (fig. 10) the frames of which are numbered to 
correspond to the numbers in the blips at the top of the pressure record 
(fig. 9) shows the position of the normal shock at various times. The 
shock moves outward between the first and second frames and then returns 
into the cowling between the third and fourth frames, where it remains 
until the fourteenth frame. Three-dimensional optical effects are 
responsible for the apparent presence of a normal shock ahead of the 
cowling in frames 4 to 13. It is evident that the flow pattern is not 
precisely symmetrical, a fact caused by a slight eccentricity in the 
location of the central body. 
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At a given value of x/L (in fig. 8) the increase of the charac- 
teristic parameters P or Q with time indicates the passage of a 
downstream or upstream compression wave, respectively. Similarly, a 
decrease in P or Q denotes an expansion wave. 

The wave pattern of the cycle is now discernible from figures 8 
to 10. The normal shock moves out because of the separation induced by 
the vortex and assumes a steady location on the spike which is dictated 
by the choking in the inlet. This motion of the shock causes an expan- 
sion wave to move downstream to reduce the mass flow. At the nozzle 
exit the downstream expansion is reflected as another expansion which 
moves upstream until it overtakes the normal shock whereupon, after a 
slight delay, this shock retreats into the cowling. The rearward motion 
of the normal shock creates a strong downstream compression wave which 
is closely followed by a much weaker expansion wave. This weaker expan- 
sion wave is probably due to the fact that when the shock moved into the 
inlet it went too far downstream. Consequently, when separation occurred 
behind it, the shock had to move up the channel somewhat and in so doing 
generated a weak expansion. In general terms, the shock overshot the 
equilibrium position when it moved back and, hence, had to adjust its 
position by moving upstream a little. It is important to note that the 
vortex from the intersection of the conical and normal shocks had to 
touch the inner wall of the cowling as the shock returned into the dif- 
fuser and yet separation did not occur to initiate the next cycle. An 
explanation for this phenomenon may be found in the following consider- 
ations. First, the normal shock is retreating so that the total pressure 
losses across it are reduced. It is then possible that inside the 
cowling the static pressure of the fluid on the outer side of the vortex 
sheet will not be required to equal its stagnation pressure; therefore, 
separation may be avoided (see reference 3). Second, the strong upstream 
expansion tends to prohibit separation by accelerating the flow and pro- 
ducing a favorable pressure gradient along the diffuser. 

The combination of a strong compression followed by a weaker expan- 
sion is reflected as a unit, with the same sign both at the nozzle exit 
and at the shock inside the cowling. Each reflection, however, lessens 
the magnitude of the waves so that, as the pressure in the diffuser is 
increased by the waves, the steps due to the waves approach zero. The 
shape of the pressure-time curve so generated has often been described 
as a "relaxation type" but is, in reality, a series of step functions 
separated by approximately constant pressure plateaus. The compression 
waves eventually force the normal shock to the cowling lip where the 
vortex attains a position from which it can start the cycle again. In 
this particular cycle, the normal shock starts to move outside the 
cowling in the beginning of the next oscillation just prior to the 
arrival of the seventh upstream reflection of the compression wave 
which by then has degenerated from an original calculated magnitude of 
4.7 pounds per square inch to only 0.1 pound per square inch. A 
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possible explanation for the origin of the second cycle between the two 
compression reflections might be that the earlier of the two waves placed 
the vortex from the normal shock in a position where a random motion 
could enable the vortex to cause separation. Another possibility might 
be a time lag required for the separation to grow to a value where it 
could choke the flow. 

The magnitude of the errors which can be introduced by the various 
assumptions may also be determined. By using the exact formula for the 
variation of the characteristic parameters across a constant-pressure 
entropy discontinuity with the lowest value of a.ja. 0oa occurring in the 

cycle, the error in the reflected value of P may be found from equa- 
tions (10), (11), and the following equations (47) and (48) (which are 
similar in form to equation (l8)): 


Qo - Q,=. - 


s 2"^e 7-1 
R 2j~ 


7 - 1 a r 


- 1, 


( 47 ) 


P 2 “ P e 



S e -S2 
— R 



( 48 ) 


This error in P is equal to 0.0019, which corresponds to a posi- 
tive pressure deviation of approximately 0.1 pound per square inch. 

AS 7-1 

« "R O'y AC y — 1 

Consequently, the assumption of « 1 and e ' % 1 + — — — 

introduces an error no larger than 0.1 pound per square inch at each 
reflection and usually is considerably less since the low temperature 
is present for only a small part of the cycle. 


The highest value of the entropy in the cycle is 0.2723; therefore, 
the maximum error involved by neglecting the passage of entropy discon- 
tinuities through the throat must not be greater than the value computed 
from equation (36), that is, Qq - Qg = -0.0059* 

This neglected increment is approximately equal to an expansion 
wave of 0.33 pound per square inch. The entering entropy, however, 
drops enough in 0.015 second after the maximum so that the net value 
of the reflected waves then would be less than 0.1 pound per square 
inch. Thus, negligence of entropy has introduced a time distortion in 
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the wave pattern by omitting an expansion wave of 0.3 pound per square 
inch followed 0.015 second later by a compression wave of 0.2 pound per 
square inch, both of which will damp out as they reflect back and forth 
in the diffuser. 


Since the strongest single compression in the cycle has a pressure 

.. « 77 . 

ratio of 

72. 

tion of isentropic waves is valid beyond question. 


= I .065 (see gage 4 of fig. 


The assumption of quasi-steady flow ahead of the cowling for pur- 
poses of entropy and mass-flow computations neglects the velocity of the 
normal shock. From the movie film it appears that the shock velocity 
is about 4 feet per second, which, computed on a normal-shock basis for 
a free-stream Mach number of I. 90 , yields a ratio of stagnation velocity 
of sound behind the shock to that of the free-stream stagnation sonic 
velocity of 1.00046 as the shock moves outward. Under these conditions, 
the tine local entropy at the lip is greater than the steady-flow value 

AS 

by an increment of = 0 . 003 , and the correct mass flow is obtained 

n 

by multiplying the steady- flow value by a factor of 0.999* The possi- 
bility of other waves between the cowling and shock must also be con- 
sidered, but, since the time interval for flow to pass from the spike 
tip to the cowling is approximately 0.0001 second and none of the 
unsteady waves has a very large pressure slope dp/dt, these waves will 
also have negligible effect. 


The influence of a variation in A e or which, of course, 

changes M e , is shown by figure 11. The A e value of 3*190 square 

inches used for the sample computations and comparison of theory and 
experiment was the channel area normal to the axis at x = 12.75 inches. 
This area is not the largest area of the diffuser, since at the end of 
the central-body support strut the area increases to 3*339 square 
inches. Thus, an effective area of 3*190 square inches presumes an 
average boundary- layer correction of 0.149 square inch in the 
remainder of the diffuser. Since there was no accurate method of meas- 
uring this correction, the influence of arbitrarily varying the effec- 

m. 

tive area to 2.99 square inches was investigated. For a — 2. = 1.000 

m 00 

and A e = 2.99 square inches, M e was found to be O.I 63 . In a similar 
manner as set forth in appendix B for A e = 3*190 square inches, the 

wave diagram and pressure-time curves were determined by using the 
identical few experimental pressures. An examination of figure 11 shows 
that the only significant difference is that the pressure in the middle 
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part of the cycle is higher for A e = 2.99 square inches than for 

ffiu 

A„ = 3.19 square inches and the same — - = 1.000. 

e m oo 

Also, the effect of changing the mass flow at start of buzz from 
^ = 1.000 to — = O.985 for an effective area A = 3*190 square 

moo ^00 

inches was determined. This computation was complicated by the fact 
that it was necessary to shift from the reflection curve at the cowling 

lin for — = 1.000 to ^ ^ = O.985 in order to complete the 

* m^ moo moo 

cycle. These reflection curves are plotted in figure 12, along with a 
line between them which was arbitrarily assumed to represent conditions 
as the shock moves from the cowling lip out to where the buzz starts 
again at a mass flow of O.985. In figure 11 the untailed diamonds 

m p 

represent the pressures if the ratio = 1.000 is assumed once the 

00 

shock has returned into the cowling until the next buzz starts. The 
tailed diamonds represent the adjustment obtained by following the 
path 1, 2, 3 (fig. 12) along the reflection curve at the cowling as the 
shock moves out. Higher pressures are evident throughout after the 
experimentally adjusted pressure portion of the cycle, with the cycle 

for — = 1.000 resulting in greater pressures than those in the cycle 
^ 00 

which allows for the mass-flow reduction from ^ = 1.000 to ^ = O.985. 

ILL qo iU CO 


A close inspection of figures 9 an( ^ H reveals the fact that the 
theoretical pressures for — — = 1.000 do not return to the original 

m oo 

trip value of 87. 0 pounds per square inch, but reach a value of only 
86.7 pounds per square inch. This discrepancy probably is due to the 
continuous slow throttling of the exit valve, since the valve-position 
indicator showed a 3*^- percent decrease in A n during the cycle from 
its value at the start of O.67 square inch. Such a throttling could 
generate upstream compression waves to increase the experimental pres- 
sure above the theoretical value computed for constant throttle area 
and Mg . 

The experimental pressure record also shows a pressure loss through 
the pipe plus an increasing hash or noise level as x/L increased. 

These results were to be expected since there were three 90° elbows in 
the standard pipe employed as a duct and in addition a chatter in the 
gate valve which was used as a throttle. The entropy increase along 
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the pipe, as well as the waves introduced by the chattering valve, of 
course, were neglected in theoretical work and are, therefore, another 
source of error. 

Approximate values of the mass flow entering and leaving the ram 
jet are presented in figure 13* The "entering" mass-flow curve is 
really an extrapolated curve determined from pressures downstream which 
gives the mass flow behind the normal shock, and, consequently, when the 
normal shock is in motion inside the cowling, an error will be present. 
This error is evident in the apparent decrease in entering mass flow 
from t = O.036 to t = 0.0375 which, in reality, is just the shock 
moving upstream inside the cowling after overshooting when it returned 
from the external spike. The mass-flow curve is therefore adjusted as 
shown by the dashed line to give the correct entering mass flow. 

An approximation for the exit mass flow was obtained by computing 
the mass flow on an effective entropy basis (square symbols in fig. 13) 
and then correcting this value when the fluid which entered the cowling 
at the start of buzz had reached the exit. Each value of the mass flow 
at the exit at time t]_ was multiplied by a correction factor based on 

the entropy existing at the cowling inlet at time t^ - At, where At 

was the time required for fluid to flow from the cowling to the exit at 
a velocity determined from the effective Mach number M e and steady 
flow. 


At = 


a,, 

M e d- 


= O.O763 sec 

0.1528 x 0.99767 x 1117.4 


( 49 ) 


Integration of the area under the mass-flow curves and division by 
the product of mass flow at start of buzz multiplied by the cyclic time 
yields the average mass-flow ratio 


J ni(t) dt 


m 

m v 


V 


(50) 


A value of 0.951 is obtained for the adjusted curve for the cowling 
inlet and 0.953 for the nozzle exit, (if a time lag At' = O.O83 second, 
which is the time required for the initial fluid to reach the exit, is 
used, a value of 0.952 results at the exit.) According to these approxi- 
mations, which include the mass adjustment shown by the shaded area in 
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figure 13, there is slightly less mass in the diffuser at the time of 
the start of the second cycle. This condition is a logical one, since 
air of higher entropy than the effective value is still flawing in the 
diffuser at time the second buzz starts. The assumptions made to com- 
pute the mass flow, however, are so rough that the accuracy of the thirc 
significant figure is very doubtful. 

The same central body and cowling were tested with the short duct 
which prescribed an over-all length of 34 inches, or ratio of length to 
maximum internal diameter of 16.5 approximating the proportions of long 
ram jets. The wave diagram and the pressure record containing both 
theoretical and experimental points are shown in figures 14 and 15 for 
the fourth oscillation after starting. This oscillation was chosen 
because of its closer similarity to the one for the longer length in 
order to facilitate comparison. Since the value of P 0 ^ f 01 * this run 

was 100.7 pounds per square inch absolute, compared to 101.0 pounds per 
square inch absolute for the longer diffuser, the value of p^ was 

reduced from 87. 0 pounds per square inch absolute to 86.6 pounds per 
square inch absolute. The theoretical computations are based on 

A_ = 3.190 square inches, — = 0.993, M e = 0.152, but the cycle is 

111 00 

mp 

not, however, adjusted for the motion of the shock from — = 1,000 
to — = 0.993 just prior to the next buzz. The value M g = 0.152 

m oo 

mu 

was chosen after an initial computation employing — = 1.000 and 

moo 

M e = 0.153 yielded results consistently low by 0.5 pound per square 

inch in the latter half of the cycle. The change in the cowling-lip 

hl mu 

parameter to reduce the mass flow from — — = 1.000 to — — = 0.993 is 

m M m oo 

less than 0° 10' (see fig. 25, reference 3) which is much smaller than 
the measurable accuracy of the parameter by the methods employed. Con- 
sequently, such an adjustment does not invalidate the solution in any 
way. 


Again, close agreement between theory and experiment is attained. 

It is evident that the compression waves computed lead the experimental 
ones by about 0.001 second. This error is due to the initial time 
location of the first adjusted compression wave being in error by that 
amount, but, since the response of the recording apparatus was flat to 
only 300 cycles per second, this error was not worth correcting. The 
valve chatter is much more pronounced in the third trace where the pres- 
sure pickup was located only 3^- inches ahead of the valve. 
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Points indicating the arrival of similar waves at gage 1 on the 
first trace of figures 9 and 15 are numbered to show the correspondence 
of the wave patterns for the long and short diffusers. The magnitude of 
the pressure variations is smaller for the short configuration because 
the reflection of the expansion wave, generated when the shock moves out 
to start the buzz, returns to hit the shock and stop its motion before 
it can move out far enough to reach a new equilibrium position on the 
spike. From figure 9 it is evident that the shock will move out for 
about 0.013 second before reaching equilibrium if its outward motion 
is due to vortex- induced separation and it is not forced out further by 
being hit by a strong compression. In contrast it takes only about 
0.005 second for a wave to traverse the short duct in both directions. 
Consequently, for the same value of A Q , all waves have less strength 

for the short diffuser than for the long one. The maximum and minimum 
pressures for the cycles computed were 87. 0 pounds per square inch and 
72.6 pounds per square inch in the long configuration contrasted to 86.6 
and 79 .O in the short one. Similarly, the minimum entering mass flow, 
referred to free-stream cowling area, was 51 percent and 79 percent for 
the long and short diffusers, respectively. 

The motion -picture film (fig. l 6 ) taken concurrently with the pres- 
sure record and numbered to correspond to the numbers on the pressure 
record shows the shock motion to be much smaller for the short pipe than 
for the long one. The shock motion also agrees with the theory in every 
respect. 

The curves of mass flow against time (fig. 17) were computed by the 
same method employed previously. The ratio of average mass flow to mass 
flow at start of buzz is 0.979 a t entrance and 0.977 a t exit. The 
entering value was not adjusted for the motion of the shock from 
m Q m-u. 

— = 1.000 to — = O .993 prior to the next cycle which would reduce 
m 00 ^ 00 

the average entering value slightly. The ratios of the average mass 

flows to the mass flow based on 100 percent capture area are 0.972 

and 0 . 970 , respectively, compared with 0 . 95 ^ and 0-953 for the longer 
diffuser. 

Since all the variations in p, a, S/R, and so forth are smaller 
for the shorter pipe, the errors introduced by the various assumptions 
(except the one made by concentrating the area changes) which simplify 
the exact one-dimensional theory are less than those for the longer pipe. 
The distortion in the wave pattern due to assuming that all the area 
changes are concentrated at the inlet and exit is increased in the 
smaller diffuser since the actual area variations occupy a larger 
percentage of the total length. (See fig. 6 .) 
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General Nature of the Buzz Cycles 

The manner in which the cyclic period is usually changed due to 
throttle adjustment is shown in figure 18, which is another part of the 
pressure record for configuration I, 0^ = 44°, L = 156 inches. At 

point I (see fig. 18 (a)) the "basic cycle" is (l) a downstream expan- 
sion, (2) an expansion reflected upstream from the exit which pulls the 
normal shock back in, (3) a downstream compression, and (4) a reflected 
upstream compression which forces the shock out to start the next cycle. 
Two smaller compressions are superimposed on the basic cycle, since the 
strong upstream compression which starts the next cycle passes through 
the contraction at the cowling entrance and reflects a weaker compression 
which shows up in points 1, 2, and 3. This is, in turn, reflected and 
re-reflected as a compression upstream from the exit, points 4 and 
downstream from the cowling contraction, points 5 "to 7> and finally 
upstream again, points 8 to 10. Further reflections are negligible. 

As the throttle is opened, the reflections from the throttle become 
smaller. This phenomenon is due to that fact that, if the conditions h, 
j, and k represent, respectively, the conditions ahead of a downstream 
wave, behind the downstream wave, and after the upstream reflection of 
the wave from a sonic valve, the ratio of the wave strengths for isen- 
tropic waves is by linearization approximately 



(51) 


As the sonic area is increased, the Mach number M^ must increase since 
the ratio of the sonic throat area to A h becomes larger and the 


reflected pressure wave becomes weaker (— approaches 1.0). 


If each reflection becomes weaker while the strength of the initial 
expansion remains approximately constant, it will take an increasing 
number of compression reflections to move the normal shock to a point 
where the vortex starts the next cycle. Figure 18(b) shows this clearly. 
Region I requires only one compression reflection; whereas regions II, 
III, IV, and V require two, three, four, and five compression reflections, 
respectively. In some cases it is evident that the cycle starts between 
reflections and, as explained earlier, starting at this time is probably 
due to the previous wave moving the normal shock to a position where 
the vortex is on the verge of causing separation so that any random 
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variation in the flow starts the cycle. Note that the basic cycle is 
approximately the same for all the frequencies and that the minor com- 
pression wave becomes unnoticeable as the valve is opened. 

The pressure record obtained when the valve was quickly opened for 
the short diffuser (fig. 19) has similar characteristics. The waves, 
however, are not as sharply defined as in the longer pipe. 

The short duct shows a marked decrease in the wave amplitude as 
the valve was opened. The magnitude of the upstream compression which 
knocks the shock out of the cowling determines the expansion strength, 
since a strong compression will knock the shock further out before it 
can be stopped by the reflected expansion. In the longer pipe, however, 
(fig. 18) the amplitude varied only slightly, since the normal shock 
moves out until it reaches an equilibrium with the sepa rat ion- choked 
cowling before any of the reflected expansion returns to alter its 
motion. 

An over-all examination of the results of the nine different 
supersonic- inlet geometries tested revealed a general behavior pattern 
for the buzzing after it was once initiated. This result was true 
although the 30° configurations were always choked at the cowling inlet 
so that the normal shock could not enter during steady flow. Also, con- 
figuration II for = 39° 50* buzzed due to the vortex from the 

separation- induced lambda shock striking the cowling. 

The pressure records for configuration III were very similar to 
those shown for configuration I 0j = 44° after the cycle in which the 

shock motion first became appreciable. This cycle was usually preceded 
by a few short oscillations of much lower magnitude. 

The pressure records for configuration II were similar to those for 
configurations I and III only after the amplitude of the buzz had become 
large enough for the basic cycle to manifest itself. With this body 
and = 43° 10’, the buzz did not start until the normal shock was 

appreciably more than halfway out on the spike and the vortex appreciably 
inside the cowling. The buzz then began in a manner best described as an 
unstable organ-pipe type of resonance. As evident from figure 20, the 
oscillations in the short diffuser increased in amplitude rapidly anrj 
then remained constant. Moving pictures showed that the buzz was not 
started by any vortex action, but rather by a boundary- layer instability, 
and that the magnitude of the oscillations built up until the normal 
shock was moving from a point one-quarter of the way down the spike to a 
point off the spike ahead of the tip. Shadowgraph excerpts (the frames 
are numbered to correspond to the numbers in fig. 20) of this configu- 
ration are shown in figure 21. The table on the following page gives 
the difference in pressures from the peak to the valley at each of the 
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three pressure pickups at times I, II, and III for three consecutive 
oscillations noted in figure 20: 



Ap for 
(psi) 


Time 

gage 1, ^ = 0.338 

gage 2, £• = 0.721 

gage 3, 2L = O.897 

1 

7.7 

10.7 

12.5 

11 

9.5 

13.1 

14.7 

hi 

11.0 

I 15.6 

18.2 


This table definitely refutes the use of a Helmholtz resonator 
analysis even for this cycle, which in the absence of starting due to 
the vortex is nearest to a sinusoidal cycle, and hence would be most 
favorable to Helmholtz theory. At the same time, it substantiates the 
plane-wave theory for the short diffuser having a fineness ratio of 
less than 17 The basis for this conclusion may be found in the primary 
assumptions employed in the Helmholtz derivation; namely, the pressure 
is uniform at any time in the large cavity, and thus the pressure fluc- 
tuations must have the same magnitude throughout the cavity. The pres- 
sure records, however, show that in the duct, which in the Helmholtz 
analysis (reference 5) is considered the resonator cavity, the magnitude 
of the fluctuations varies by about 50 percent in the same cycle. ur- 
thermore, although the oscillation is unstable and is increasing rapidly, 

the pressure variation near the valve (%■ = O.897) during one cycle is 

still larger than the pressure variation further upstream ^ = 0.338) 

during the next cycle. In other words, the diffuser behaves similarly 
to an organ pipe with relative nodes in velocity and anti-nodes m pres- 
sure at the exit for this cycle. Additional verification was obtained 
from the longer version of the same supersonic-diffuser configuration 
where a typical oscillation generated a pressure variation of 1.2 pounds 

per square inch at ^ = 0.074; whereas at £ = 0.705 the variation was 
4.5 pounds per square inch. 

Figure 20 also is an exception to the general trend of increasing 
frequency with decreasing valve area since the periodic time increases 
from 0.0149 second to 0.0155 second during resonance. The reason for 
this result is found in the "inertia" effects of the external motion of 
the shock. From figures 14 and 16 it is obvious that the shock does 


NACA TN 3695 


33 


not start- to return to the cowling as soon as struck by the upstream 
expansion wave (t = 0.0114 second, frame 5) but delays before starting 
inward (t = 0.0152 second, frame 7). As the amplitude of the shock 
motion increases, more of the expansion will have to strike it to 
initiate its movement toward the cowling. Thus, for figure 20, the 
periodic time, comprised of the lag time and basic cycle time (that is, 
that time required for an expansion and compression to traverse the duct 
in both directions), increases as the amplitude and lag increase. 

It is then obvious that, for any finite-amplitude oscillation gov- 
erned by reflections from the exit, the periodic time must be greater 
than the basic cycle interval. The minimum period for the configura- 
tions tested was 0.0148 second, with a basic cycle time of 0.0101 sec- 
ond, when waves reflected from the exit controlled the buzz. Configu- 
ration II, Q z = 36° 50’, however, oscillated at a very high frequency 

with a period of approximately 0.002 second when the exit area was 
reduced appreciably. This oscillation may have been governed by the 
supersonic diffuser entirely and perhaps resulted from the periodic 
boundary- layer variations passing through the first minimum. 


A Possible Application to Simulate Constant-Area Heat Addition 


By Cold Flow With Area Contraction 

Once the quasi-one-dimensional nature of the wave phenomena has 
been established, the possibility arises of employing area changes to 
simulate heat addition. An examination of equations (l) and (2) shows 
that the parameters P and Q are functions only of the changes in A 
or S/R. It is possible therefore to get the same variation for one 
(not both simultaneously) of these parameters due to a change in S/R 
while A is held constant that can be obtained by holding S/R constant 
and varying A. A generalization of this principle is used to produce 
theoretically a similar wave pattern for the case of constant-area heat 
addition and isentropic area contraction for configuration pairs like 
those shown in figure 22. 

Consider the simplified combustion process in a semi-infinite pipe 
diagrammed in figure 22(a) where the combustion is assumed to occur in 
a discontinuous flame front. If conditions at the flame front may be 
assumed to reach a steady or equilibrium flow after each transmission 
and reflection of a pressure wave, the equations of continuity, energy, 
momentum, and state may be incorporated in a trial-and-error process to 
determine the waves transmitted and reflected for each incident wave. 

(The exact characteristic equations cannot be used since the derivatives 
become infinite at the discontinuity. ) Of course, an assumption regarding 
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the amount of energy added per unit mass flow must be made for the 
energy equation. This assumption is a very difficult assumption to 
make since the waves change the Reynolds number, velocity, and air- fuel 
ratio at the combustion zone. 

Two very simple assumptions were used for the two illustrative com- 
putations; namely, (l) the ratio of total temperature across the flame 
front was constant, and (2) the energy added per unit mass was constant. 
For an initial Mach number of 0.200 ahead of the flame and 0.400 behind 
it, these assumptions prescribed a constant stagnation-temperature ratio 
of 3.0482 for the first case and an increment in stagnation temperature 
of 2.0482 times the initial, undisturbed, stagnation temperature for the 
second case. A downstream-moving expansion was then assumed to strike 
the flame and theoretical computations made for the wave pattern due to 
the reflections and re- reflections from the flame discontinuity and 
nozzle. Under the limitation of a choked exit nozzle the waves rapidly 
damped out. The general form of the wave diagram is shown by the solid 
lines of figure 23 and the corresponding pressure-time curve at x = 

(located ahead of the flame by an amount equal to one-half the distance 
from the flame to the nozzle) for three different initial expansion 
strengths is plotted in figure 24 . 

Computations were then made for a discontinuity in area in place of 
the flame front with the requirements that: (l) undisturbed conditions 

ahead of the discontinuity be identical to those for the heat addition, 
and (2) the medium- strength expansion have the same first upstream 
reflection for both the area change and heat addition. By trial and 
error the area ratio corresponding to the constant-stagnation-temperature 
ratio was found to be 0. 5129 and that for the constant increment in stag- 
nation temperature to be O.38OI. The length of the system after the 
discontinuity in area change was determined so that an infinitesimal 
disturbance would require the same amount of time to travel from the 
discontinuity to the exit nozzle and back to the discontinuity for both 
the heat addition and contraction cases. In other words, if e a denotes 

the distance from the flame front to the nozzle and e ~ the distance 
from the area change to the nozzle (see fig. 22), 
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The rest of the wave pattern until damping was completed was then 
computed under the same restrictions as imposed on the heat-addition 
solution. Good agreement is shown in figure 2k for the pairs of 
pressure-time curves for the corresponding conditions of (l) total 
temperature ratio equal to 3*0482 and area ratio of 0.5129^ or (2) total 
temperature increment of 2.0482 and area ratio of O.38OI. For each of 
the three initial expansions, all the curves except that for the con- 
stant at o will damp to the same final pressure (which is, of course, 

different for each expansion) and return to the original Mach number 
of 0.200 ahead of the discontinuity. Thus, for either of the heat- 
addition assumptions made, a corresponding area change could be employed. 

When the amplitude of the buzz in an actual ram jet is small, the 
variation of Reynolds number, velocity, and air-fuel ratio is also small. 
In a ram jet, however, there is also the added complication of coupling 
between the fuel- injection nozzles and the combustion zone whereby a 
wave originating from the combustion zone travels upstream with a 
velocity a - u to the fuel-injection system and changes the air-fuel 
ratio of the fluid present at the fuel nozzles. This fluid then travels 
downstream with a velocity u to the combustion zone where the new air- 
fuel mixture burns to generate a new wave. This phenomenon cannot be 
duplicated by waves whose speed is strictly limited to the local sonic 
velocity superimposed on the local fluid velocity. If the reflections 
from the shock at the inlet, from the nozzle exit, and from the waves 
striking the flame front, however, have a stronger influence on the 
pressure pattern than the waves generated by the change in air-fuel 
ratio, there is the possibility that the combustion in a particular 
buzzing ram jet may follow approximately a law of energy release which 
can be simulated by an area change. In that case the actual ram jet 
(fig. 22) with a finite combustion zone might be tested for its behavior 
during buzz by the cold- flow test of the configuration (fig. 22) in which 
the combustion chamber is replaced by a shorter section with a gradual 
area reduction. 


CONCLUDING RIMARKS 


Experimental instantaneous pressure records of simulated ram-jet 
models without heat addition showed that these models buzzed in a manner 
governed internally by quasi-one-dimensional (plane-wave) unsteady-flow 
theory and not in accord with an analysis of the problem by the laws of 
Helmholtz resonator theory. A theory applicable to the buzzing problem 
was obtained through modification of the one -dimensional, unsteady-flow 
theory by introducing specified assumptions which greatly facilitated 
the computational efforts without impairing the usefulness of the theory 
for this flow. Because of the difficulty in computing the three- 
dimensional, unsteady, viscous flow when the shock was in motion ahead 
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of "the cowling and in computing when the shock was in motion far from 
equilibrium with the waves moving upstream from the nozzle, it was 
necessary to use a few experimental points at the start of the cycle as 
a basis on which to compute the remaining pressure points of the cycle. 
Theoretical computations of the pressure-time curves, which were made 
on this basis for one supersonic-diffuser configuration with two differ- 
ent subsonic-duct lengths, agreed closely with the experimental records. 
The entrance and exit mass-flow - time curves, as well as the wave dia- 
grams, were also computed for these configurations by the unsteady-flow 
theory. 

The wave pattern causing the buzz was found to be a basic cycle 
which consisted of a downstream expansion reflected at the exit as an 
upstream expansion and a downstream compression reflected as an upstream 
compression plus a variable number of compression-expansion couplets 
following the basic cycle. The number of compression-expansion couplets 
depended on the valve-exit area and decreased as this area decreased; 
therefore, the frequency generally increased as the throttle was closed. 

The simulation of combustion in' a ram jet by an area contraction 
for buzzing tests was discussed as a possible extension of the results . 


Langley Aeronautical Laboratory , 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , January 29, 1952- 
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APPENDIX A 

DERIVATION OF CHARACTERISTIC EQUATIONS 


The general form of the quasi -one-dimensional continuity equation 
'is expressed as follows.’ 


A(pA) + J^-(puA) = 0 (54) 


where the area A is considered to he a function of hoth time and posi- 
tion. With this equation, cases of a boundary layer growing with time 
as well as continued throttling may be treated. Equation (5*0 may be 
written in another form as 


d logp A d log_ A d log„ . p d log„ p 

-S — + u — + — + u - h — 

dt dx dt dx dx 


0 


( 55 ) 


The following two identities 


= 7 RT 


( 56 ) 


n S = _d_(heat added) 
R _ RT 


( 57 ) 


are now introduced. The quantity a in equation ( 56 ) is the speed of 
an isentropic wave of infinitesimal strength, or sound wave, which may 
move through a fluid of varying entropy. The M heat added 11 of equa- 
tion ( 57 ) includes the heat due to convection, dissipation, and sources. 
When the energy equation 


d(heat added) = dT - ^ dp 


( 58 ) 


and the equation of state 


p = pRT 


( 59 ) 
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are combined with equations (56) and (57), the differential equations 
between the properties of state and the nondimens ional entropy are 
determined 


d I = TT d l0g e a 


d logg p 


(60) 


d 




log e 


a 


d log e p 


(61) 


Although these relationships were derived for a particular elemental 
volume of fluid, they also are true for all the various parts of the 
fluid and may be used when going from one part of the fluid to another. 
This procedure is justified because the entropy is a property of state 

so that S(x, t) = s[p(x,t), p(x,t), a(x,t)] if all the fluid was origi- 

nally in the same state. Consequently, the following partial derivatives 
may be obtained: 


N x S 

d logg p _ 2C V d log e a 0 R 

ax r ax ax 


(62) 


6 logg P = 20y. a logg a - 6 | 

at r at at 


(63) 


6 loge 

ax 


n a s 

2Cp a logg a R 


r a x 


ax 


(64) 


a logg 
at 


p _ 


2C p a logg a ^ R 


r at 


at 


(65) 


The momentum equation 
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becomes equation ( 67 ) upon insertion of equations ( 5 6 ) and (59) 


Su + u du + a£ a log e , P = 0 (67) 

St Sx 7 Sx 


Substitution of equation (64) into equation ( 67 ) and simplifying yields 

S 


du du 2Cy da _ a^ a R . 

dt + U dx R dx 7 dx 


( 68 ) 


C P 

since 7 = 7T-. 

W 

Similarly, substitution of equations (62) and ( 63 ) into equa- 
tion (55) results in: 

8 log- A 8 log_ A 2C„ 8 2C V is- a r An 

a — + au + -^2a-a^ + — ^u2a-au,--£i+a2ii=0 

dt dx R dt dt R dx dx dx 


( 69 ) 


Addition of equations ( 68 ) and ( 69 ) and then rearrangement yields 

( 2 C v \ , , d /2C V \ 

— - a + u + (u + a)li- a + u = 

V R / dx\ R / 


_d_ 

dt 


/d logg A + d log e A j + 7 


R 


dt 


dx 


<bt 


+ u 


8 |\ 


d | d | 

— 

+ a 

— - + (u + a)— S. 
St ox 

dx / 

7 



( 70 ) 


Similarly, subtraction of equation (68) from equation ( 69 ) and rearrange- 
ment gives 
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iL a - u) + (a - aO- & 
Stl R / Sx\ R 


a - u = 


IdJfBeJ. + u + L^l a(il ♦ u it) + & 

\ st Sx / y \st Sx / 7 


The following identities are used: 


(l) Convective derivative + u 

v ' Dt St ox 


i-S + (u 

St 


a)|i 

Sx_ 

( 71 ) 


(72) 


dx 

(2) Total derivative along a characteristic of slope — = 


u 1 a, 


A = <L + ( u + a)A 

5t St Sx 

^ .. ?r., r „ 

(3) P s 

(4) Q EE 


(73) 


7 - 1 a o 00 

p| 

0 ! 

8 1 

11 

R a O 00 a Ooo 

2 a 

u _ 

20 v a u 


7 " 1 a o 00 a o ro R a o 


00 ^00 


Insertion of these terms in equations (70) and (71) gives the character- 
istic equations in their final form 


s s 

§p a D lp 6e A , 7 - 1 a D R + 1 a S J| 
8t ” a 0oo Dt 7 a^Dt 7\6t 


s s 

5Q _ a D lQ ge A , 7 - 1 a D R , 1 a 5 R 
St "a 0oo Dt 7 a 0oo Pt 7 a o 00 8t 
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APPENDIX B 

COMPUTATION FOR CONFIGURATION I, = 44°, SUPERSONIC DIFFUSER 
WITH OVERALL RAM-JET LENGTH OF 1 56 INCHES 


The first step in the computation for configuration I, = 44°, 

supersonic diffuser with overall ram-jet length of 156 inches is the 
solution of the following equation (which is the same as equation (A3) 
of the text) 


M e ( 1 + L ^ M e 2 


1 

2 




7+1 

2(7-1) 


For the effective area A e , the value of the area at x = 12.75 inches 

is chosen; this area is equal to the area of the annulus formed by the 
inside of the pipe, which forms the principal part of the diffuser, and 
by the supporting strut for the central body which extends to 
x ** 15.5 inches (fig. 5) 


A e = j(2.062 2 - 0.436 2 ) = 3.190 


(72) 


Also, since to the order of accuracy of the no-flow shadowgraph 
measurements (fig. 25 of ref. 3), the buzz should start at a mass -flow 

ratio ^ = 1.000 because the conical shock falls on the cowling lip. 

When these values are used with 7 = 1.400 plus the measured values 

p 

of Aj_ = ^(1.2) square inches, p e = 87. 0 pounds per square inch 
absolute (fig. 9), p n = 101 pounds per square inch absolute (recorded 

4-' 00 

during test ) , and the nozzle Mach number = l-90> equation (^*3) can 

be solved for M e = 0.1528. The effective entropy -1- (based on ^ = 0 

is determined from the assumption of steady flow just prior to the start 
of buzz; thus, since the above data yield an effective stagnation pressure 
in the duct at start of buzz of 88.4 pounds per square inch absolute 


£*b _ s_ 

R R 


88.4 


^ - -log e ^ = 0 * U32 


101 


(73) 
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Similarly, 

1 


2 

f£_ = (1 + 7 ~ 1 x 0 . 1528 2 ) = 0.99767 

aooo V 2 / 

(74) 

Ub - 0.1528 X 0.99767 = 0.15244 

a °oo 

(75) 

which yield starting values of the parameters P and Q 


P b = 5.1^08 

(76) 

Q b = 4.8359 

(77) 

The reflection curve (solid line in fig. 25) for the entrance is 
computed as set forth in the theoretical section by adjusting the quasi- 
steady values of P 2 and Qg b 7 a factor which is a function of the 

c 

difference — - 0.1332. Reflections at the nozzle exit are specified 
R 

5 - M 

by the relationship ^3 = r + y[ ^3 0. 9^-069? 3 (dashed line in fig. 25) 

The final curve (fig. 26) needed for the completion of the 

P ( a/a °co V/ 

is the pressure - sonic-speed curve Q ^ 99767 / ' see 

wave diagram 
equation 21 ) . 


Experimental values, designated by the flagged symbols on the pres- 
sure record (fig. 9 ) and wave diagram (fig. 8 ), are used to determine the 
downstream-moving waves when the normal shock is in rapid motion or ahead 
of the cowling. A linear interpolation of the variation of the param- 
eters P and Q, is assumed for all waves. 

The remainder of the wave diagram is constructed by the two reflec- 
tion curves and then the pressures found on the third curve from the local 

, .. a _ P + Q 

sonic velocity Jq — • 


Approximate values of the entropy at the entrance may be found from 
the quasi- steady flow to the inlet when the shock is ahead of the cowling. 
Consider the conditions existing (fig. 8 ) at t = 0.0193 the entrance 
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where P = 5.OO6I and Q = 4.8359 so that = O.O85I, r 2 - = 0.98420 

°oo °oo 

(in the imaginary lamina) and p = 79*10. The quasi- steady- flow energy 

equation gives the actual = 0.999275 and, since at a free-stream 

a °oo 


Mach number of I.90, 


15.07 for p 

U 00 


101.0, and 

’ a n 


0. 76205, 



0.999275 

0.98420 


!°g e 


79.10 

15.07 


0.2390 


(78) 


The linearized approximation might he used instead to obtain the value 


s 2 ~ S oo _ s 2 ~ S e + s e ~ S oo _ 2y Aa 
R R R 7 - 1 a 

- 7 0-9993 - 0.9 8 4 2 + 0<1332 . 
0.9993 


+ 0.1332 

0.1058 + 0.1332 = 0.2390 (79) 


The ratio of the instantaneous entering mass flow to the mass flow at 
start of buzz is then, from equation (29): 


^2 = /5.0061 + 4.8359x5.0061 - 4.8359 A _ 0. 1058 \ 
m b \5.l4o8 + 4.8359/5.1408 - 4.8359 \ 3.5 / 


= O.506 


( 80 ) 
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(a) Typical ram jet. 





(b) Simplified ram jet assumed for one -dimensional, unsteady-flow 

analysis. 

Figure ,1.- Sketch of typical ram jet and simplified one -dimensional 

counterpart. 


Cowling 
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Figure 2.- Method employing quasi-steady 
cowling inlet to obtain value of P^' 


flow values ^ 2 , Q 2 ) a_ t 

reflected from shock by 


Q 
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(a) Variation of parameters across an entropy discontinuity 
where pressure and velocity are constant. 

Figure 3 .- Diagrams showing how characteristic parameters P and Q 
may he referred to effective entropy by insertion of imaginary 
lamina of fluid in the flow field. 
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(b) Insertion of imaginary lamina of fluid (c) with entropy into 

flow (d) of different entropy returns parameters to value existing 
originally in state (3) to permit calculation of pressure and 
velocity. 


Figure 3-- Concluded. 


Cowling 
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Figure 4.- Wave generated by passage of entropy discontinuity through 

exit nozzle. 


Nozzle 


VJ1 

o 



Figure 5.- General configuration of models tested. 

in inches. 


All dimensions are 
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(b) Medium diffuser; L = 156 inches. 


O 



(c) Long diffuser; L = 390 inches. 


Gate valve 


Location of pressure pickup , n (not necessarily same for various 
diffusers) 

P Location of 90° elbows 


0 .1 .2 .3 .4 .5 .6 .7 

Position along ramjet axis, x/L 

Figure 6.- Schematic diagram of ram-jet models tested showing: 

(1) relative length of central body to overall length, and 

(2) locations of pressure pick-ups and elbows. 



008 0 


VJ1 

ro 


0 ° 36 ' 



Figure 7.- Central body of configuration I. (Dimensions are in inches). 
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Figure 8 .- Wave diagram for configuration I. L = 156 inches. 

Flagged symbols denote points computed from experimental pressure 
record. 


Pressure, psia 





Figure 9 .- Theoretical and experimental pressure-time curves for con 1 
figuration I. 9 Z = 44°; L = 156 inches. Flagged symbols denote 
experimental points on which theoretical computations are based. 
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Figure 10 .- Low-speed shadowgraphs during buzz of configuration I. 

= lpip 0 ; L = 156 inches. 



Pressure , psia 
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(a) Curves for gage 1 at ^- = 0.07^. 

Figure 11.- Theoretical pressure-time curves showing effect of variation 
in — and A e . Flagged symbols denote common pressures taken from 

m oo 

experimental record. 
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(b) Curves for gage 4 at — = 0.558* 


Figure 11.- Concluded 
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Figure 12.- Curves for determination of parameters reflected from inlet. 

g 

Effective entropy -§■ = 0.133^* 

K 


Instantaneous mass flow 

Mass flow at start of cycle 



Figure 13 . - Curves of instantaneous mass flow through inlet cowling and 
exit nozzle for configuration I. 9^ = 1<4°; L = 156 inches; 


moo 


is assumed to he equal to 1 . 000 . 
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Figure 14. - Wave diagram for configuration I. 0 1 = L = 3^ inches. 

Flagged symbols denote points computed from experimental pressure 
record. 


Pressure , psia 



Figure 15 *- Theoretical and experimental pressure-time curves for con- 
figuration I. = Ml- 0 ; L = 3 ^* inches. Flagged symbols denote o\ 

experimental points on which theoretical computations are based. 


NACA TN 3695 


n n n n n in 


62 


NACA TN 3695 


I c c c 

• c c c 

• C- C- C" 
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• c- 

• c- c- c 

L-72716 

Figure 16.- High-speed shadovgraphs during buzz of configuration I. 

Q z = 14.14°; L = 3 1 *- inches. 



Instantaneous mass flow 
Mass flow at start of cycle 



.004 .008 .012 .016 .020 .024 .028 .032 .036 .040 .044 .048 .052 

Time, sec 


Figure 17-- Curves of instantaneous mass flow through inlet cowling and 
exit nozzle for configuration I. 0 1 = 44°; L = 34 inches; 

mK 

jjj— is assumed to be equal to 0.993- 
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(a) Enlargement of region I. 

Figure 18.- Pressure record obtained as exit throttle is rapidly opened 
for configuration I. = 44°; L = 156 inches. 
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(b) Complete record. 
Figure 18.- Concluded. 
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Figure 19-- Pressure record obtained as exit throttle is rapidly opened 
for configuration I. = M)-°; L = Jh inches. 



Figure 20.- Pressure record showing unstable oscillation as buzzing 
begins for configuration II. = ^3° 10'; L = 3^ inches. 
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Figure 21.- High-speed shadowgraphs during buzz of configuration II. 

9 j = k3° 10 * ; L = 34 inches . 
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(a) Discontinuous heat addition in semi-infinite pipe. 



h *£ >1 


(b) Discontinuous area change in semi-infinite pipe. 



(c) Ram jet with finite combustion- chamber length. 



(d) Ram jet with gradual contraction of rear section. 

Figure 22.- Configuration pairs for possible simulation of heat 

additions. 
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Characteristics for heat addition 

Characteristics for area change 


Figure 23 


Simplified wave diagram for heat addition and corresponding 
area contraction in semi- infinite pipe. 


Nozzle for 
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Figure 2k. - Comparison of pressure-time curves for discontinuous heat 
addition and area contraction in semi-infinite pipe. 
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Downstream characteristic parameter, Pj or P ‘ 3 

Figure 25-- Curves for determination of parameters reflected from 
cowling inlet and nozzle for configuration I. 0^ = 44 , 

L = 156 inches. 
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Figure 26.- Curve of local static pressure and local sonic speed for 
configuration I. L = 156 inches. 
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